Critical temperature of a Bose gas in an optical lattice 
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We present theory for the critical temperature of a Bose gas in a combined harmonic lattice potential based on 
a mean-field description of the system. We develop practical expressions for the ideal-gas critical temperature, 
and corrections due to interactions, the finite-size effect, and the occupation of excited bands. We compare our 
expressions to numerical calculations and find excellent agreement over a wide parameter regime. 
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Ultra-cold atoms in optical lattices have emerged as a flexi- 
ble system for studying many -body physics (QlSlltl- Despite 
immense interest in this system, our understanding of the fun- 
damental phenomenon of Bose-Einstein condensation, in par- 
ticular the critical temperature, Tc, is rudimentary compared 
to the case of the harmonically trapped gas. Much of the diffi- 
culty in making predictions for the lattice system arises from 
the complex spectrum ifl IH 01 in the combined harmonic lat- 
tice (CHL) potential used in experiments. In contrast, the 
power-law density of states for systems in purely harmonic 
traps has led to simple expressions for Tc and the finite- 
size and (mean-field) interaction Igt] corrections. Although 
in combination these corrections are at the ~10% level for the 
harmonically trapped gas, beautiful experiments by the Or- 
say group llOil have been able to exclude ideal-gas behavior 
by more than two standard deviations, and find quantitative 
agreement with mean-field theory. The prospects for simi- 
lar thermodynamic studies in optical lattice experiments has 
taken a leap forward with the recent demonstration of a prac- 
tical scheme for measuring temperature in this system ifTTll 
(see also fll ). 



Here, we first derive an expression for the ideal-gas crit- 
ical temperature, T^", of the experimentally relevant CHL 
system by introducing a shape approximation to the lattice 
band structure. By explicitly including the physics of the 
finite width of the ground band, our result demonstrates re- 
markably good agreement with the exact (ideal) result 
We then derive expressions for the critical temperature cor- 
rections due to finite size, excited band, and mean-field in- 
teraction effects. We compare our expressions to full self- 
consistent three-dimensional (3D) numerical calculations us- 
ing the actual band structure and demonstrate their excellent 
agreement over a broad range of experimentally relevant pa- 
rameters. We observe that the interaction correction is most 
significant, heavily suppressing Tc. We find that the validity 
range of the quadratic shape approximation is complementary 
to the effective-mass approximation, so that simple descrip- 
tions extend over a wide range. All of our expressions de- 
pend only on experimental parameters and quantities that can 
be simply evaluated from the single-particle spectrum of the 
translationally-invariant lattice (TIL). 

The general CHL potential we consider is Viatt (r) + Vtr (r), 
for the case of 1 < d < 3 spatial dimensions. The lattice po- 
tential is given by Viatt(r) = ^ {'^'''j /'^j)' where 



Vj is the lattice depth and Oj is the lattice site spacing, in di- 
rection 7. The harmonic trap is Vtr (r) = 5"^ ^j''^' 
where m is the atomic mass and ujj is the harmonic trap fre- 
quency. We also consider the TIL case, where Vtr (r) = 0. We 
note that while all the numerical results we present are for the 
cubic lattice case (all aj ~ a, Vj ~ V), our expressions also 
apply to non-cubic lattices. 

The Hartree-Fock description of bosons in an optical lattice 
at T > Tc gives the per-site density of atoms in band b ifisll 
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We have used the local density approximation, valid for kT ^ 
to*, where Co* = u* /d,uj* = luj ^m/m* and m* is the 

ground-band effective mass along direction j (evaluated at 
quasi-momentum). Here, gb{K) is the density of states of the 
corresponding TIL (i.e. usual Bloch spectrum) and the inter- 
action coefficient Ubb' is proportional to the s-wave scatter- 
ing length fls and the integral of the product of Wannier state 
densities for bands b and b' fl4[l (also see lITsll ). We have also 
introduced (3=1 /kT and the chemical potential, /i. Equation 
([T]i, derived from an extended Bose-Hubbard Hamiltonian, 
becomes ambiguous in the shallow lattice limit (V^ < 2Ej^) 
where the Wannier states become delocalized. We implement 
a continuation of the theory into this regime (see also Ill4fl ). 
but note that very shallow lattices are best examined without 
making a band decomposition |ITS|. We also note that mean- 
field theory is invalid for interacting systems in the deep lattice 
limit where correlations are important The breakdown 

of mean-field theory has been studied for a one dimensional 
optical lattice at T = fl^l- For experiments with ®^Rb, ([B 
is a good description for V < 13Ej^. 

The total number of thermal atoms in band b is Nt = 
a^** f dr n(,(r). For a wide parameter regime of interest, the 
critical point occurs when the occupation of excited bands is 
negligible. Initially, we consider the ground band {b = 0) 
in isolation and assume it contains all atoms in the system 
(i.e. iVo = N). We will later account for the influence of 
excited bands. 

The bandwidth is an important parameter, which for band 
b, we calculate as 
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FIG. 1: Exact (solid) and quadratic shape approximation (dash) 3D 
cubic TIL ground-band density of states for (a)V = and (b) V 
5Er, where V = V^"^^ Er = h'^/Sma^ and a = Hi 



l/d 



where K™™ is the minimum energy of the band. This defi- 
nition for Wb provides a useful characterization of the band- 
width in shallow lattices (e.g. see Fig.[T]i, approaches the tight- 
binding expression in deep lattices (see Fig.|2ja)), and extends 
our treatment to non-cubic lattices. 

Ideal-gas critical temperature. — For sufficiently deep lat- 
tices, the TIL ground bandwidth {W = Wq) becomes neg- 
ligible compared to other relevant energy scales. Then, the 
exact structure of the ground band becomes unimportant and 
band shape approximations can be made. As a first approxi- 
mation, we set go (A') = S{K — W/2) /a'^. Neglecting interac- 
tions, ([T]) gives nQ{r) = cxp[/314r(r) +w/2] — l} ^, 
where we have defined the thermal bandwidth w = W/kT, 



and z 



Then iVo = NX 
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where C,^{z) 



^^-^ z^'^ /n°' is the Bose function, d is the number of har- 
monically trapped dimensions (i.e. d = for the TIL), Ns is 
the number of sites for the TIL and Ng = {2'KkT /muj'^a^Y/'^ 
is a measure of the thermally accessible number of sites for 



the CHL, with w ee ]\^ uj^'^. In the limit w ^ d, for 
d < 3, or for the TIL, there is then no condensation, and for 
the 3D CHL, the critical temperature is 
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with ({a) = Ca(l)i as in |6tl (also see |[l9ll ). A better approxi- 
mation to go{K) is to use the quadratic shape approximation, 
go{K) = 6K{W - K)/W^a'^ for < < 1^ and zero 
otherwise |20tl, as shown in Fig. [T] Using ([1]), this gives 
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For the 3D TIL this resuh gives T^" w WiN/N, + ^)/3k UM- 
For the 3D CHL, the series expansion of Q, convergent for 

w < 27r is 1I2T 
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Solving this implicit Taylor series, we find the ideal-gas crit- 
ical temperature (with = W/kT^', the notation for the 
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FIG. 2: (a) Width of the 3D cubic TIL ground band, using ^ (solid), 
maximum less minimum energies (dash) and the tight-binding ap- 
proximation Ql, I23I1 (dot), (b) Spread function S{w) using l ll2t 
(solid) and the quadratic shape approximation with il It (dash). 

thermal bandwidth, uses the subscripts and superscripts from 
the corresponding temperature) 
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Interaction shift. — With harmonic confinement, the lead- 
ing order effect of interactions on Tc is at the mean-field level 
in 3D: repulsive interactions reduce the density at trap center 
and decrease Tc. Our approach, similar to involves a Tay- 
lor series expansion about the critical point. However, in the 
lattice, the id eal g as has divergent density at r = as w ^ 0, 
from (HI and 1I22I1 . The effect of a non-zero interaction param- 
eter, U = Uqo, is significant at removing this near-divergent 
behavior. Thus, making the Taylor expansion at T° would sig- 
nificantly overstate the interaction effect, particularly for deep 
lattices. We therefore use the first terms in a Taylor series 
about the interacting critical temperature, Tc, (the quantity we 
aim to determine) and obtain, to first order in 1^1 and U 

STc = Tc- « -i^i^^No + UduNo)/dTNo, (9) 

where UdjjNo ~ —2U /dr no(r)9^no(r)/a'^ and /i 
2f/no(0), with the derivatives {dx = ^) and no(0) evaluated 
at /.t = 0, [/ = 0, and T = T^. We obtain 

STc « -2C/[1 - Siwc)]noWd^No/dTNo, (10) 

where Wc = W/ kTc and we have defined the spread function 

' Jdr no(r)9,,no(r) 
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fj.=a 



Using the approximation go{K) = S{K — W/2) / a^, we find 
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For the quadratic shape approximation, as w ^ 00, S{'w) 

+ 't)"^/VC(2)C(f ) « 0.325, and we have 
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evaluated S{w) numerically as shown in Fig. |2|b), which 
shows that 0.24 < S{w) < 0.325 and, as u; ^ 0, S{w) -> 
0.24. We note that for the (no lattice) 3D harmonic trap case, f.^- 
5 = 0.281 H. 

We evaluate the derivatives required by ( fTOl i analytically, 
then expand using ll22ll to find (to first order in U) 
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where the last result uses the approximation S{wc) 
giving a simple equation to be solved for Tc- 

Finite-size effect. — Accounting for the non-zero energy 
of the ideal ground state (and hence the saturated ideal-gas 
chemical potential) due to confinement, eq = 3huj* /2, gives 
rise to the so called finite-size shift in the condensation tem- 
perature . To quantify this effect for the 3D CHL we obtain 
6T'^' « -eod^No/drNo yielding 
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where w° = W/kT^ 
binding effective mass, we have uj 
and eo 



For a cubic lattice, using the tight- 

UJTTy/W/12Eli m 



^/SW/ErttHuj/A where Cu = '^j^j/d- In this 
limit, the highest order finite-size shift in (fl4] i depends on W 
only through T° 
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We compare the finite-size effect using the quadratic shape 
approximation to the full numerical calculation (from ([T) with 
/i = eo, limiting the integral to A'+Vtr(r) > eo) in Fig-Ha.b). 

Excited bands. — Over a wide regime of experimental in- 
terest, the excited bands are weakly occupied at the critical 
point (i.e. kTc < K™™ for b > 0). It is then sufficient to ac- 
count for the small occupation using a simple band shape for 
excited bands, gb{K) = l/Wha'^ for < if - ifj"'" < Wb 
and zero otherwise, and neglecting excited-band interactions. 
We find the total excited band occupation (which acts to re- 
duce the ground-band occupation from N) at temperature T 



b>0 



W,/feTj^ (17) 



where we have used the approximation C^(e ^) « e ^ for 
x > 1. If iVcx(T'c) < N (noting Tc is often much smaller 
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FIG. 3: (Color) 3D cubic CHL results. (a,b) Comparison of ideal 
finite-size effect between quadratic shape approximation (solid) and 
numerical results (ground band in green and all bands in blue), 
effective-mass l il8b (dash-dot) and ui — > limit l il6t (dash). (c,d) 
Comparison of quadratic shape approximation (solid) to numerical 
results (dot), using a — 425 nm, as — 5.77 nm and ignoring the 
finite-size effect. Both ideal (ground band in green, all bands in blue) 
and interacting results (ground band in red, all bands in magenta) are 
shown. For these parameters, the Mott insulator transition for unit 
occupancy is at V Ri 13Er. (a,c) cj = 27r x 32 Hz, = 10^. (b,d) 
= 27r X 64 Hz, iV = 10*'. 

than T^ (fT3])). we may safely neglect excited bands. If excited 
band contributions are important, it is necessary to use a con- 
sistent procedure to determine Tc. The approach we use is: 

1. Calculate an initial estimate T° using ([H]) if w° ^ 2n (oth- 
erwise, solve (|5])). 

2. Calculate the ground-band interaction shift, STc, using ( flj] ). 

3. Adjust the ground-band number No = N- Nc^{T° + STc). 

4. Recalculate r° using No- 

Finally, steps 3 and 4 are iterated to find T° if necessary. The 
critical temperature allowing for interactions in the ground 
band, with an adjustment for excited bands is then + STc. 

Results. — We have verified the usefulness of our expres- 
sions by comparison with full self-consistent numerical calcu- 
lations using the actual band structure, ([TJ, over a wide range 
of experimentally relevant parameter regimes. The ideal-gas 
results in Fig. [3ta,b) reveal that the finite-size correction is 
typically only a few percent of T°, and is usefully bounded by 
the limiting expressions ( fT6l ) and ( fTSl l. Fig. [3|c,d) examines 
the more significant interaction and excited band effects. The 
system in Fig. Oc) is typical of current experiments and re- 
veals the large suppression of Tc from r° due to interactions, 
and negligible contribution from excited bands. The system 
considered in Fig. Od), with larger atom number and tighter 
harmonic confinement, has higher critical temperature giving 
appreciable excited band effects |25}|, especially in the shallow 
lattice regime {V < 4:Ef>). We observe, from the full range of 
results in Fig. [3] that our expressions are in good agreement 
with the numerical calculations. 

Effective mass. — In the low temperature regime {kTc <C 
W) the quadratic shape approximation is less effective as 
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FIG. 4: Comparison of actual (dot) to effective-mass (dasli) and 
quadratic sliape approximation (solid) critical temperature for an 
ideal Bose gas in a 3D cubic CHL with (a) a; = 27r x 16 Hz, 

N = 10^ (b) Lj = 27r X 64Hz, = 10^ 

the detailed features of go{K) at low K become important. 
However, for an ideal gas in this regime we can employ 
an effective mass approximation. For the TIL case, this 
predicts no condensation for d = 1,2, and yields Tj^ = 

{2TTh'^/km*) [N/Nsa^Ci5/2)]^^'^ in 3D IH. For the CHL, 
there is no condensation for d ~ 1 and for d = 2, 3, we have 
= hw* [N/C{d)f''^ /k H. We compare the effective- 
mass predictions to our full numerical calculation for a 3D 
CHL in Fig. SJa), showing good agreement in the low tem- 
perature regime, although the effective-mass approximation 
is usually very poor except for very shallow lattices (e.g. 
Fig. Sfb)). The finite-size effect for the 3D CHL, using the 
methodology of 11], is 



TP 



C(2) 



oj* 2C(3)2/3 



7V^i/3 _o.73^tLiv-i/3, (18) 



as shown on Fig.[3la,b) (for a cubic lattice Q* /oj* ~ uj/lo). 

Conclusions. — We derived practical expressions for the 
ideal and mean-field interacting critical temperature in an op- 



tical lattice by using a shape approximation for the density 
of states. We have given primary consideration to the exper- 
imentally relevant case of the 3D combined harmonic lattice 
potential, but have also given results for the translationally- 
invariant lattice. For the ideal case we have shown that our 
validity range (low to moderate values of w) is complemen- 
tary to that of the effective-mass approximation (large values 
of w). By using the average energy of the band to calculate the 
width in (|2]i, we have extended the validity range into shallow 
lattices. 

We have compared all of our results to full self-consistent 
numerical calculations with the actual band structure and have 
generally found excellent agreement over a wide parameter 
regime. We have derived a formula for the finite-size effect in 
the lattice, and the important effect of mean-field interactions 
in the ground band. We have also demonstrated a procedure 
to include the influence of excited bands. 

Our results indicate that, for the typical parameters of cur- 
rent optical lattice experiments, interaction shifts in Tc are 
large compared to the harmonically trapped case (where the 
interaction shift is usually only a small fraction of T° lloi fioll ). 
and that excited band effects are small. Our main results, con- 
tained in (|7]i, ( fT3] l and ( fT4l ). are easy to evaluate, depending 
only on experimental parameters and the width of the ground 
band. Our theory is applicable to non-cubic lattices with ap- 
propriate identification of the bandwidth parameter (|2|. With 
the recent progress in measuring temperature in lattices we 
believe our results will be directly comparable to experiments 
in the near future. 
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